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High-pT particle spectra in p + p (p + p), p + A and A + B collisions are calculated within a 
QCD parton model in which intrinsic transverse momentum, its broadening due to initial multiple 
parton scattering, and jet quenching due to parton energy loss inside a dense medium are included 
phenomenologically. The intrinsic kr and its broadening in p + A and A + B collisions due to 
initial multiple parton scattering are found to be very important at low energies ( y/s < 50 GeV). 
Comparisons with S + S, S + Au and Pb + Pb data with different centrality cuts show that the 
differential cross sections of large transverse momentum pion production (pr > 1 GeV/c) in A + B 
collisions scale very well with the number of binary nucleon-nucleon collisions (modulo effects of 
multiple initial scattering) . This indicates that semi-hard parton scattering is the dominant particle 
production mechanism underlying the hadron spectra at moderate pr ~ 1 GeV/c. However, there 
is no evidence of jet quenching or parton energy loss. Within the parton model, one can exclude 
an effective parton energy loss dE q /dx > 0.01 GeV/fm and a mean free path \ q < 7 fm from the 
experimental data of A + B collisions at the SPS energies. Predictions for high pr particle spectra 
in p + A and A + A collisions with and without jet quenching at the RHIC energy are also given. 
Uncertainties due to initial multiple scattering and nuclear shadowing of parton distributions are 
also discussed. 

PACS numbers: 25.75.-q, 12.38.Mh, 13.87.-a, 24.85.+p 



I. INTRODUCTION 

Large- Et partons or jets are good probes of the dense 
matter formed in ultra-relativistic heavy-ion collisions 
[1,2], since they are produced in the earliest stage of 
heavy-ion collisions and their production rates are cal- 
culable in perturbative QCD. If a dense partonic matter 
is formed during the initial stage of a heavy-ion collision 
with a large volume and a long life time (relative to the 
confinement scale 1/Aqcd), the produced large Et par- 
ton will interact with this dense medium and, according 
to many recent theoretical studies [3-5], will lose its en- 
ergy via induced radiation. The energy loss is shown to 
depend on the parton density of the medium. Therefore, 
the study of parton energy loss can shed light on the prop- 
erties of the dense matter in the early stage of heavy-ion 
collisions. It is also a crucial test whether there is any 
thermalization going on in the initial stage of heavy-ion 
collisions. 

Even though one cannot measure directly the energy 
loss suffered by a high energy parton propagating through 
a dense medium, its effective fragmentation functions 
must change and leading particles must be suppressed 
due to the parton energy loss. This can be measured 
directly in deeply inelastic e + A collisions (for jets go- 
ing through a normal nuclear matter) [6,7] and direct- 
photon-tagged jets in high-energy heavy-ion collisions 
[7,8]. Since the large pr single-inclusive particle spec- 
tra in nuclear collisions are a direct consequence of jet 



fragmentation as in hadron-hadron and hadron-nucleus 
collisions, they are also shown [9] to be sensitive to parton 
energy loss or modification of jet fragmentation functions 
inside a dense medium. 

Because of the extremely small cross sections of large 
Pt particle production at the CERN SPS energy (a/s « 
20 GeV), measurements of particle spectra in heavy- 
ion collisions at moderately large pt have not become 
available before recent experiments [10,11]. In a re- 
cent paper [12], the author used a QCD parton model 
calculations to analyze the large pr spectra and found 
that there is no evidence of parton energy loss in high- 
energy heavy- ion collisions. From the analysis of the 
spectra, one can limit the effective energy loss to less 
than dE/dx < O.OlGeVVfm. Since experimental data 
have shown evidences of an interacting hadronic matter 
in heavy-ion collisions at the CERN SPS energy, one can 
also conclude that the hadronic matter dose not cause 
apparent jet energy loss. 

In this paper, we will conduct a systematic study of 
high pt hadron spectra in pp, pA and AB collisions, 
since the reliability of the determination of a small ef- 
fective parton energy loss in high-energy heavy-ion col- 
lisions crucially depends on the precision to which we 
understand the spectra without parton energy loss. We 
will include both the intrinsic fc^ m p + p collisions and 
Ut broadening due to multiple initial-state scattering in 
p + A and A + A collisions in a phenomenological man- 
ner, which we find necessary to describe large pt particle 
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production at energies below y/s < 50 GeV. We will then 
compare the calculations with the experimental data for 
S + S, S + Au and Pb + Pb collisions at the CERN SPS 
with different centrality cuts and verify the scaling be- 
havior characteristic of hard processes. Using the same 
parton model we will then predict the high pt spectra 
in pA collisions at the RHIC energy and then calculate 
same spectra in AA collisions with and without effects of 
parton energy loss. The goal of this paper is to provide 
a practical and phenomenological description of the high 
Pt hadron spectra in p + p and p + A collisions. This 
will provide a baseline prediction for high-energy A + A 
collisions against which we can study the effects of the 
parton energy loss caused by the dense medium. 



II. P + P COLLISIONS: INITIAL K T 



Large px particle production in high-energy hadron- 
hadron collisions has been shown to be a good test of the 
QCD parton model [13]. Assuming partons inside the col- 
liding hadrons only have longitudinal momenta, the large 
Pt particle production cross section can be calculated 
as a convolution of elementary parton-parton scattering 
cross sections, parton distributions inside the hadron and 
parton fragmentation functions. The factorization theo- 
rem ensures that the parton distributions and fragmenta- 
tion functions are universal and can be measured in other 
hard processes, e.g., e + e~ annihilation, deeply inelastic 
ep collisions. With high precision data, parton distribu- 
tions and fragmentation functions have been parameter- 
ized including their QCD evolution with the momentum 
scale Q of the hard processes [14,15] which in turn can be 
used to calculate particle production in many other hard 
processes. This parton model with collinear approxima- 
tion has been rather successful so far in describing large 
Pt particle production in high-energy pp collisions [16] 
with > 50 GeV. 

The above described collinear parton model fails to 
account the experimental data on angular correlation of 
produced heavy quarks and the total transverse momen- 
tum distribution of the heavy quark pairs [18] or the 
Drell-Yan lepton pairs [17]. In terms of parton mod- 
els one naturally expects that partons inside a hadron 
carry at least an average intrinsic transverse momentum 
of about a few hundred MeV, reflecting the hadron size 
via the uncertainty principle. Furthermore, higher order 
pQCD processes like initial state radiation or 2 — > 3 sub- 
processes with additional radiated gluons can also give 
rise to large initial transverse momentum for the collid- 
ing partons. However, even within next leading order 
(NLO) pQCD collinear parton model, experimental data 
of Dell-Yan lepton pairs production [17], heavy quark 
pair production [19], and direct photon production [20] 
all indicate that an average initial parton transverse mo- 
mentum of about 1 GeV is needed to describe the data. 
Such a large value is an indication of its perturbative 



nature. Theoretical techniques have been developed to 
include all high order contributions [21-23], though it 
is still problematic to differentiate what is true intrin- 
sic and what is pQCD generated initial transverse mo- 
mentum. One can re-sum all the radiative corrections 
into a "Sudakov" form factor up to a scale Qq and de- 
fine (non-perturbative) contributions below this scale as 
"true" intrinsic, even though such an "intrinsic" average 
transverse momentum still depend on colliding energy 
y/s and the momentum scale Q of the hard processes. 
Such schemes have been developed for processes, such as 
Drell-Yan [21-23] and heavy quark production [24] where 
two scales are involved, with only minimum phenomeno- 
logical success in limited kinematic region. The problem 
in processes such as direct photon and single inclusive 
hadron production, where there is only a single scale, is 
even more difficult [25]. 

In this paper, we will adopt a more phenomenolog- 
ical approach by introducing initial parton transverse 
momentum with a Gaussian form of distribution as has 
been done in Drell-Yan [26], heavy quark [19], direct pho- 
ton [20] and single inclusive hadron production [20,27,28]. 
Since it is difficult to extend the NLO calculation to p+A 
and A + A collisions where multiple collisions occur and 
our main purpose is to provide an effective description 
of single inclusive hadron spectra not only in p + p but 
also in p + A and A + A collisions, we will use lowest or- 
der (LO) parton model with a simple scheme of including 
the intrinsic transverse momentum with a Gaussian dis- 
tribution. The inclusion of initial transverse momentum 
can be considered as a parameterization of both higher 
order and nonperturbative corrections. The predictive 
power of this model then lies in the energy and flavor de- 
pendence of the hadron spectra. Similar comments can 
be made for our modeling of multiple parton scattering 
effect in p + A and A + A collisions. 

One can define parton distributions in both fractional 
longitudinal momentum x and initial transverse momen- 
tum kT in a factorized form, 

dxd 2 k T gN(kT,Q 2 )fa/N(x,Q 2 ), (1) 

where gN(kT,Q 2 ) is the transverse momentum distribu- 
tion and f a /N{ x , Q 2 ) ar e the normal parton distribution 
functions for which we will use the MRS D-' parame- 
terization [14]. The inclusive particle production cross 
section in pp collisions will then be given by [13] 
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where D^ c (z c , Q 2 ) is the fragmentation function of par- 
ton c into hadron h as parameterized in [15,29] from e + e~ 
data, z c is the momentum fraction of a parton jet car- 
ried by a produced hadron. The K w 2 (unless otherwise 
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specified) factor is used to account for higher order QCD 
corrections to the jet production cross section [30]. 

We define the momentum fraction x in terms of the 
light-cone variables x a = (E a + p\\ a )/\fs for partons in 
the forward beam direction and Xb = (Eb — p\\ b )/Vs in 
the backward beam direction. The four-vector momenta 
for the colliding partons are then p a = (E a , kx a ,P\\a)> 
p b = (E b ,\c T b,P\\b) witn 

E a = (x a ^+^)/2, 

Xay/S 

p lla = (x a VS-^=)/2, 
E b = (x by ^+-^=)/2, 

Xby/S 

P\\ b = (-x b V~s+-^=)/2. (3) 



ffJV (fc T ,Q 2 ) = _l_e- fe -/< fe -K (5) 



Since our initial kx includes both the intrinsic and pQCD 
radiation-generated transverse momentum, the variance 
in the Gaussian distribution (k T )^ should depend on the 
momentum scale Q of the hard processes. As pointed 
out by Owen and Kimel [28], in addition to the leading 
log contribution from initial-state radiation that is in- 
cluded in the Q 2 dependence of the parton distributions, 
the recoil effect of large-angle gluon emissions can lead 
to Q-dependent (h^N in this approach of QCD parton 
model. In the following we shall consider both a constant 
(k T )N and one that depends on the momentum scale of 
the hard processes. For the Q-dependent case, we choose 
the following form [28] in our model 



In principle, one should also include the transverse mo- 
mentum smearing from the jet fragmentation. We ne- 
glect this in our calculation of particle spectra in the 
central rapidity region and consider it been effectively 
included in the calculation by adjusting the initial kx dis- 
tribution. In this case, particles are produced in the same 
direction of the fragmenting jet. The Mandelstam vari- 
ables for the elementary parton-parton scattering pro- 
cesses are then, 

k 2 k 2 

s = x a Xb\fs + Ta Th - 2k Ta k Tb cos(0 a - <p b ), 

%a%b\/ S 

i = -— (x a sf S e-y + -^ey - 2k Ta cos <j> a ) , 

Zc V X W S J 

u = - V — (x b ^se y + -^=e- y -2k Tb cos <h) , (4) 

Zc V X aV S J 

where px and y are the transverse momentum and rapid- 
ity of the produced particle, cos</> a = k/r a ■ Pr/kxaPr, 
and cos^b = kx b ■ px/kx b px- The momentum fraction 
of the produced hadron z c is then given by the identity 
for massless two-body scattering, s + t + u = 0. 

In Eq. (2), one should also restrict the initial trans- 
verse momentum kx a < x a\/s and kxb < x b ^fs such that 
the partons' longitudinal momenta in Eq. (3) have the 
same signs as their parent hadrons. In addition, one of 
the Mandelstam variables can approach to zero if the 
initial kx is too large and then the parton cross section 
could diverge. To avoid these problems, we introduce a 
regulator p 2 in the denominators of the parton-parton 
cross sections da/di(ab — ► cd). We choose /j, — 0.8 GeV 
in our following calculations. The resultant spectra are 
sensitive to the choice of p only at pr around p where 
we believe that QCD parton model calculation becomes 
unreliable. 

We assume the initial fcy distribution (fer) to have 
a Gaussian form, 



{k 2 T ) N (Q 2 ) = 1.2(GeV 2 ) + 0.2a s (Q 2 )Q 2 - (6) 



The form of the Q-dependence and the parameters are 
chosen to reproduce the experimental data. Following 
the same approach as in Refs. [27,28], we choose Q 2 to 
be Q 2 = 2siu/(s 2 + P + u 2 ). 

Shown in Figs. 1-4 are our calculated spectra for 
charged pions as compared to the experimental data 
[31,32] for p + p collisions at E lah =200, 300, 400, 800 
GeV. As one can see from the figures that QCD par- 
ton model calculations with the initial kx smearing (solid 
lines) as given in Eq. 6 fit the experimental data very well 
over all energy range. However, without the initial kx 
smearing (dot-dashed lines for tt~) the calculations sig- 
nificantly underestimate the experimental data, as much 
as a factor of 20 at £q a b = 200 GeV. This is because 
the QCD spectra are very steep at low energies and even 
a small amount of initial kx could make a big increase 
to the final spectra. As the energy increases, the QCD 
spectra become flatter and small amount of initial kx 
does not change the spectra much, as we can already see 
by comparing the spectra for £i a b = 800 GeV (Fig. 4) to 
that of 200 GeV (Fig. 1). This is further demonstrated 
in Fig. 5 where we compare the LO calculation to exper- 
imental data in high-energy pp collisions. Here we have 
used a K = 1.5 factor. One can sec that at higher collider 
energies, the initial kx does not make much difference to 
the spectra at high transverse momentum. In the case of 
no intrinsic kx, we didn't use the regulator p in our cal- 
culation. Instead, we used the px of particle as a cut-off 
of the phase integral in Eq. (2). This is why the spec- 
tra (dot-dashed lines) increase faster at smaller px than 
the ones with initial kx smearing in which a regulator of 
p = 0.8 GeV is used. 
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Slab = 



2 3 4 5 6 7 8 
p T (GeV/c) 

1. Single- inclusive pion spectra in p + p collisions at 
200 GeV. The solid lines are QCD parton model cal- 
culations with Q-dependent intrinsic kx and the dot-dashed 
line is without. The dashed line is for a constant intrinsic 
(&§-)jv = 1.5 GeV 2 . Experimental data are from Ref. [31]. 
The inserted figure shows the corresponding tv~ /tt + ratio. 




3 4 5 6 
p T (GeV/c) 

FIG. 2. Single-inclusive pion spectra in p + p collisions at 
Slab = 300 GeV. The solid lines are QCD parton model cal- 
culations with Q-dependent intrinsic kr and the dot-dashed 
line is without. Experimental data are from Ref. [31]. The 
inserted figure shows the corresponding tv~ /tt + ratio. 



In Fig. 1, we also show as a dashed line the calculated 
spectrum of ir~ with a constant average initial trans- 
verse momentum, (fcf.) at = 1.5 GeV 2 . At large pt such a 
constant intrinsic (fcf^jv underestimates the experimen- 
tal data. This demonstrates phenomenologically why the 
Q dependence of {k 2 -.) n is needed to fit the experiment- 
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2 3 4 5 6 
p T (GeV/c) 

FIG. 3. The same as Fig. 2, except at E Xah = 400 GeV 




2 4 6 8 10 
p T (GeV/c) 
FIG. 4. The same as Fig. 2, except at £"i a b 
the experimental data are from Ref. [32] 
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800 GeV and 



al data in particular at large values of xt — 2p T /^/s for 
the choice of Q 2 we used. One can also use an alter- 
native choice for Q 2 = Pj?* = (pt/z c ) 2 . In this case 
we found that one docs not need Q-depcndcnt (h^N 
to fit the data in pp collisions. However, one needs to 
use an energy-dependent initial parton transverse mo- 
mentum distribution. We also checked that as energy 
increases such a Q-dependent {k 2 -) n is not needed any- 
more to fit the experimental data. However, the available 
experimental data at collider energies are only limited to 
finite pt range where Xt = 2pt/\/s is not so large as 
compared to the low energy data. In the following we 
will use the Q-dependent initial momentum distribution 
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FIG. 5. Single-inclusive spectra of charged hadrons in p + p 
collisions at y/s = 200, 900, 1800 GeV. The solid lines are 
QCD parton model calculations with intrinsic kr and the 
dot-dashed lines are without. Experimental data are from 
Refs. [33,34]. 
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FIG. 6. Single-inclusive spectra of charged kaons in p + p 
collisions at i?i a b = 200 GeV. The lines are QCD parton model 
calculations with intrinsic kr and data are from Ref. [31]. 

for our study of high px spectra in p + A and A + A 
collisions. 

In the inserted boxes in Figs. 1-4, we also plot ir~ /n + 
ratio as a function of px- At higher px, particle pro- 
duction is more dominated by the leading hadrons from 
valence quark scattering. Since there are more up-quarks 
than down-quarks in p + p system, one should expect the 
ratio to become smaller than 1 and decrease with px- The 
QCD parton model calculations describe this isospin de- 



pendence of the spectra very well. For other flavors of 
hadrons, we show in Figs. 5-7 the calculated kaon spec- 
tra at different energies. We see that the agreement with 
experimental data is very good, except at large px at low 
energies which may be improved by using better fragmen- 
tation functions for kaons. One can also notice that K~ 
spectra at large pr are about a factor of 10 smaller than 
K + . This is because the content of strange quarks in a 
nucleon is much smaller than up (down) quarks which 
are responsible for leading K + (K°) hadron production. 
We did not calculate and compare the spectra of pro- 
duced protons and anti-protons because there is no very 
accurate parameterization of the corresponding fragmen- 
tation functions. 
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FIG. 7. The same as Fig. 6, except at E ub = 300 GeV 
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FIG. 8. The same as Fig. 6, except at E ub = 400 GeV 
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III. P + A COLLISIONS: K T BROADENING 

Since single-inclusive particle spectra at high pt in 
p + A collisions have been shown, both experimentally 
[31] and theoretically [35], to be sensitive to multiple 
initial-state scattering, or Cronin effect, it is important 
that we take into account this effect here in order to have 
a quantitative study of the change of particle spectra at 
high px in nucleus-nucleus collisions. One can study the 
Cronin effect in a model of multiple parton scattering 
[35]. In such a model, one uses Glauber multiple scat- 
tering formula to treat parton scattering between beam 
and target partons which is essential to incorporate the 
interference effect. For example, when contribution from 
double scattering is considered, one must also include 
the absorptive part of the single scattering which has a 
negative contribution. This absorptive part then cancels 
part of the contribution from double scattering. Conse- 
quently, the enhancement of particle spectra because of 
double scattering decreases with px in the form of 1/p? 
and in general with [36] . One can also find that the 
dominant contribution in double scattering comes from 
the case where one of the scatterings is soft and large 
part of the final px of the produced jet comes from just 
one hard scattering. At much lower pr the absorptive 
correction is much larger than the double scattering con- 
tribution and the spectra there is even suppressed. So the 
integrated cross section or average particle multiplicity at 
moderate px are not affected by multiple parton scatter- 
ing even though the spectra are modified, as observed 
in the DY case [37]. This will provide us a justification 
for our following phenomenological treatment of multiple 
parton scattering in p + A and A + A collisions. 

In this paper we assume that the inclusive differen- 
tial cross section for large px particle production is still 
given by a single hard parton-parton scattering. How- 
ever, due to multiple parton scattering prior to the hard 
processes, we consider the initial transverse momentum 
kx of the beam partons is broadened. Assuming that 
each scattering provide a kx kick which also has a Gaus- 
sian distribution, wc can in effect just change the width 
of the initial kx distribution. Then the single inclusive 
particle cross section in minimum-biased p + A collisions 
is, 



dyd 2 p T 



K 



abed 



2 bt A (b) J dx a 



dxbd k a xd kbx 



x 9A(k aT , Q 2 ,b)f a/p (x a ,Q 2 ) 

x g P (k bT , Q 2 )h/A(xb, Q 2 , b) 

Dl /c {z c ,Q 2 )do. 

(ab — > cd), 



irz c 



dt 



(7) 



where tA(b) is the nuclear thickness function normalized 
to J d 2 bt A (b) = A. We will use Wood-Saxon form of 
nuclear distribution for tA(b) throughout this paper un- 
less specified otherwise. The parton distribution per nu- 



cleon inside the nucleus (with atomic mass number A and 
charge number Z) at an impact parameter b, 



fa/A(x,Q 2 ,b) = S a/A (x,b) 



A 



fa/p(x,Q 2 ) 



(l-^)f a/n (x,Q 2 ) 



(8) 



is assumed to be factorizable into the parton distribu- 
tion in a nuclcon .f a /N{ x iQ 2 ) an d the nuclear modifica- 
tion factor S a /A(x, b) which we take the parameterization 
used in HIJING [38] for now. The initial parton trans- 
verse momentum distribution inside a projectile nucleon 
going through the target nucleon at an impact parameter 
b is then, 



gA(k T ,Q 2 ) 



1 



n(k T ) A 
with a broadened variance 



(9) 



(k 2 T ) A (Q 2 ) = (k 2 T ) N (Q 2 ) + 5 2 (Q 2 ){v A {b) - 1). (10) 

The broadening is assumed to be proportional to the 
number of scattering v A (b) the projectile suffers inside 
the nucleus. For the purpose of considering the impact- 
parameter dependence of the kx broadening, we will sim- 
ply assume a hard sphere nuclear distribution. Therefore, 



v A {b) = a NN t A (b) = a NN 



3A 
2irR 2 A 



l-b 2 /R% (11) 



where R A = 1.12A 1 / 3 fm and aNN is the inelastic 
nucleon-nucleon cross section. 

We also assume that kx broadening during each 
nucleon-nucleon collision S 2 also depends on the hard 
momentum scale Q — J^f* , the transverse momentum 
of the produced parton jet. Such a dependence is easy to 
understand by comparing the role of in the multiple 
parton scattering case to that in the initial and final state 
radiation associated with a hard jet production. For ex- 
ample, unless the multiple scatterings are far separated 
in space-time such that the propagating parton becomes 
on-shcll, the interaction cross of the propagating parton 
will depend on its virtuality which in return depends on 
the scale of the hard scattering. In Rcf. [39], the kr 
broadening in a multiple parton scattering scenario is re- 
lated to the ratio of the leading-twist parton distributions 
and higher-twist multiple parton distributions which ap- 
parently depends on the momentum scale of the hard 
parton scattering. In this paper, we will use the follow- 
ing scale-dependent kr broadening per nucleon-nucleon 
collision, 



5 2 {Q 2 ) = 0.225 



ln 2 (Q/GeV) 
1 + ln(Q/GeV) 



GeV 2 /c 2 



(12) 



Such a functional form is chosen to best fit the experi- 
mental data in p + A collisions. For Q = 2 ~ 3 GeV, 



G 



S 2 = 0.064 ~ 0.129 GeV 2 /c 2 , which is consistent with 
the value obtained from the analysis of kx broadening 
for J/\I> production inp + A collisions [40,41]. We should 
point out that this value can be different from what one 
gets from the analysis of kx broadening of Dell-Yan lep- 
ton pairs in p + A collisions, where mainly quarks and 
anti-quarks are involved in the hard processes and there 
is no collision effect after the Dell-Yan production point. 




10 2.5 

p T (GeV/c) 

FIG. 9. Ratios of charged pions spectra in p+W over p+Be 
each normalized by the atomic number of the target nucleus. 
The lines are the parton model calculation with kr broaden- 
ing due to multiple parton scattering. The dot-dashed line 
in the third panel is the calculated result with a constant 
kr broadening S 2 = 0.1 GeV 2 . Experimental data are from 
Refs. [31,42]. 

Using the above assumptions, we calculate the large px 
particle spectra in p + A collisions and compare to the 
spectra in p+p collisions. As we show in the last Section, 
particle spectra has a strong isospin dependence in the 
parton model. In order to minimize this known isospin 
dependence of the spectra in our study of nuclear depen- 
dence of the spectra, we compare the spectra of p + A 
for heavy target with that for a very light nuclear target. 
Shown in Fig. 9 are our calculated ratios of charged pion 
spectra (solid lines for ir + and dot-dashed lines for ir~) 
in p + W over that of p + Be each normalized by the 
atomic number of the target nucleus. The experimental 
data are from Ref. [31,42]. If there was no nuclear depen- 
dence due to multiple scattering, the ratios would have 
a flat value of 1, modulo the residue isospin dependence 
because of the small isospin asymmetry of the target nu- 
cleus. The difference between our calculated ratios for 
ir + and 7r~ gives the order of this isospin asymmetry ef- 
fect. One can get rid of this effect by using the ratios of 
7r + + 7r~ spectra. As shown in the figure, our model can 
roughly describe the general feature of the nuclear depen- 
dence of the spectra at large px due to multiple parton 



scattering. The ratios should become smaller than 1 at 
very small pt because of the absorptive processes as we 
have mentioned. But here our perturbative calculation 
will eventually breaks down because of the small mo- 
mentum scale. At larger pt, the spectra are enhanced 
because of multiple parton scattering. As px increases 
further, the ratios decrease again and saturate at about 
1. The decrease follows the form of l/p T consistent with 
the general features of high twist processes. Since the 
transverse momentum broadening due to multiple par- 
ton scattering is finite, its effect will eventually become 
smaller and disappear. Therefore, the px location of the 
maximum enhancement can give us the scale of average 
transverse momentum broadening. We also show in the 
third panel as dot-dashed line our calculation with a con- 
stant kx broadening, 5 2 — 0.1 GeV 2 . While the result 
has the general feature of the experimental data, the fit 
is not as good as the scale-dependent kx broadening. 
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FIG. 10. Ratios of charged kaons spectra in p + W over 
p + Be each normalized by the atomic number of the target 
nucleus. The solid (dot-dashed)lines are the parton model 
calculation for K + (K~) with kr broadening due to multiple 
parton scattering. Experimental data are from Refs. [31,42]. 

At low energies, E\ a b = 200 GeV for example, the 
full structure of multiple scattering cannot be revealed 
because of the finite phase space constrained by the ki- 
netic limit. Furthermore, the ratio will increase at large 
xx = 2px/\/s because of the so-called EMC effect [43] 
on nuclear structure function at large x which is incorpo- 
rated into the parameterization of the nuclear modifica- 
tion factor for parton distributions inside a nucleus. One 
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word of caution should also be given here. Our calcu- 
lation might not be reliable anymore near the boundary 
of the available phase space (pr ~ y/s/2) where overall 
energy and momentum conservation is very important 
which is not imposed in QCD parton model. 

We have also calculated the kaon spectra in pA colli- 
sions. Show in Fig. 10 are the ratios of the kaon spectra 
mp + W over p + Be collisions (normalized by the atomic 
number). The agreement with the experimental data is 
quite good for K + . However, for K~ there is quite sizable 
discrepancy between the calculation and the data. This 
is might be due to the flavor dependence of the kr broad- 
ening suffered by the quarks, ft is quite possible that u 
quarks have larger scattering cross sections as they prop- 
agate through the nucleus and therefore experience larger 
kx broadening than s quarks. Such effects will result in 
the observed difference in nuclear enhancement between 
K~ and K+. 



IV. A + A COLLISIONS: JET QUENCHING ? 

It is straightforward to incorporate the initial kr 
broadening due to multiple parton scattering in A + A 
collisions. Since partons from projectile and target 
beam both suffer multiple scattering before the hard pro- 
cess, their initial transverse momentum distributions are 
broadened as given by Eq. (9). In addition, the pro- 
duced parton jets have to propagate through a dense 
medium and interact with other produced particles in 
the medium. Theoretical studies [3-5] show that a fast 
parton propagating inside a dense partonic medium will 
suffer considerable amount of energy loss. Since large 
Pt particles are produced through the fragmentation of 
parton jets with large transverse energy, parton energy 
loss will definitely lead to the suppression of large pr 
particles. One can describe this suppression phenomcno- 
logically via effective fragmentation functions which are 
modified from their original forms in the vacuum [7-9]. 
In this approach one can then calculate the single inclu- 
sive particle spectra at large px in A + B collisions as, 



d<J AB 

dycPpT 



K^J d 2 bd 2 rt A {r)t B (\h-v\) 

abed 

x J dx a dx d 2 k a xd 2 k T 

x 9A{kaT,Q 2 ,r)g B (k b T,Q 2 , |b - r|) 
x .f a /A{xa,Q 2 ,r).f b / B (xb,Q 2 , |b - r|; 
D h/c (z c ,Q 2 ,AL) da 



dt 



(ab cd), (13) 



where D h / C (z c , Q 2 , AL) is the modified effective fragmen- 
tation function for produced parton c which has to travel 
an average distance AL inside a dense medium. Wc will 
not elaborate on the modeling of the modified fragmen- 
tation functions [7-9] here, except pointing out that it 



depends on two parameters: the energy loss per scat- 
tering e c and the mean free path A c for a propagating 
parton c. The energy loss per unit length of distance is 
then dE c /dx = e c /X c . We also assume that a gluon's 
mean free path is half of a quark and then the energy 
loss dE/dx is twice that of a quark. In principle, the 
energy loss dE/dx should depend on local parton den- 
sity or temperature, the parton's initial energy, the total 
distance and whether there is expansion [44] . These pos- 
sible features will influence the final hadron spectra and 
their phenomenological consequences have been studied 
in detail in a previous publication [9]. In this paper, we 
will simply assume a constant energy loss and study the 
average effect of parton energy loss in dense matter. 

After integration over the impact parameter space, 
Eq. (13) will give us the inclusive hadron spectra for 
minimum-biased events of A+B collisions. If one neglects 
nuclear effects (parton shadowing, k T broadening and jet 
quenching), the resultant spectra should be exactly pro- 
portional to AB which is the averaged number of binary 
nucleon-nucleon collisions. Very often experimentalists 
also measure inclusive hadron spectra for certain classes 
of events with different centrality cuts (total transverse 
energy Et or charged multiplicity in the central region) . 
Since one cannot directly measure the impact parameters 
in each class of centrality, a theoretical model of corre- 
lation between impact parameter b and the total trans- 
verse energy Et (or charged multiplicity) has to be in- 
troduced in order to calculate the inclusive cross section 
for events with different centrality cuts. In this paper, we 
will use the correlation function P AB (ET,b) [normalized 
to / dE T P AB {E T ,b) = 1] introduced in Ref. [45] which 
assumes a wounded nucleon model for transverse energy 
production. The differential cross section 



da ab 
dE T 



= J d 2 b[l 



<T A B W] p AB{ETjb) (14) 



has been shown to reproduce the experimental data 
of NA35 and NA49 [46] very well, where T AB (b) = 
J d 2 rt A (r)t B (b — r) is the overlap function for A + B 
collisions at impact parameter b. We refer readers to 
Ref. [45] for the details of this model. 

After incorporating the correlation function between 
impact parameter and transverse energy in Eq. (13), the 
inclusive cross section of large pt hadron production in 
events with centrality cut E T e (£™ in , EJp**) can be 
shown to be roughly proportional to 



(^binary) = tfb / dE T P AB (E T , b)T AB (b) (15) 

which is just the average number of binary nucleon- 
nucleon collisions. For minimum-biased events, this num- 
ber is just AB. Since the actual da AB / dEx distribution 
depends on each experiment's coverage of phase space 
and the detector's calibration, we will choose the values 
of E™™ and _E™ ax so that the fraction of the integrated 
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cross section in Eq. (14) within the Et range matches 
the experimental value of a given centrality cut. 
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FIG. 11. Single-inclusive spectra of 7r in S + S, S + Au and 
Pb + Pb collisions (both minimum-biased and central events) 
at the CERN SPS energies. The solid lines are QCD parton 
model calculations with kr broadening due to initial multiple 
parton scattering and the dashed lines are without. Experi- 
mental data are from Refs. [10,11]. 

Shown in Fig. 11 are the calculated inclusive spectra 
for produced ttq in S + S, S + Au and Pb + Pb collisions, 
both minimum-biased and central events. The QCD par- 
ton model calculations with the kx broadening due to ini- 
tial multiple scattering (solid lines) agree with the exper- 
imental data (WA80 and WA98) [10,11] well very at p T 
above 1 GeV/c. No parton energy loss has been assumed 
in the calculations. The dashed lines are the spectra in 
pp collisions at the same energy multiplied by the nuclear 
geometrical factor as given in Eq. (15). The difference 
between the solid and dashed lines is simply caused by 
effects of kx broadening and nuclear modification of par- 
ton distributions inside nuclei. These effects are similar 
as in p + A collisions and are more important in colli- 
sions of heavier nuclei. Without these nuclear effects the 
high pt hadron spectra in A + B collisions are exactly 
proportional to the average number of binary collisions 
as shown by the dashed lines. This is a common charac- 
teristic of hard processes in p + A and A + B collisions. 
Because of absorptive processes, low px particle produc- 
tion, which can be considered as coherent over the di- 
mension of nuclear size, has much weaker A-dependence. 
In the wounded-nucleon model, soft particle production 
cross section is proportional to the average number of 
wounded nucleons which is much smaller than the num- 
ber of binary collisions. By studying the transition of 
the scaling property of the hadron spectra from low to 
high pt values, one can then determine at what pr value 
the underlying mechanism of hadron production become 



dominated by hard processes. This will be the subject of 
another study. 




p T (GeV/c) 

FIG. 12. Parton model calculations of single-inclusive spec- 
tra of 7T° with different values of parton energy loss dE q /dx 
and mean free path X q in central Pb + Pb collisions at 
E\ a b = 158 AGeV/. The Experimental data are from 
Ref. [11]. 

In the parton model calculations presented so far, ef- 
fects of parton energy loss have not been considered yet. 
If there is parton energy loss and the radiated gluons be- 
come incoherent from the leading parton, the resultant 
leading hadron spectra at large px from the parton frag- 
mentation should be suppressed as compared to p+p and 
p + A collisions. As we have shown, however, the parton 
model calculations without parton energy loss fit the ex- 
perimental data very well. To find out the consequences 
of an effective parton energy loss in the hadron spec- 
tra and how the experimental data constrain its value at 
the CERN SPS energy, we show in Fig. 12 as the dot- 
dashed line the calculated ttq spectra in central Pb + Pb 
collisions with dE/dx — 0.01 GeV/fm. In the calcula- 
tion (dot-dashed line) we have assumed the energy loss 
per scattering e q — 0.02 GeV and the mean free path 
X q = 2 fm (dE q /dx — e q /X q ). In our model [7,8], the 
modification of the fragmentation function is sensitive to 
both parameters. The mean free path is a measure of 
the strength of the interaction between the leading par- 
ton and medium while energy loss dE/dx reflects the de- 
gree of attenuation induced by the interaction. Shown as 
the dashed line in Fig. 12 is a calculation with the same 
dEq/dx — 0.01 GeV/fm but with a mean free path X q = 7 
fm which is about the average total length a parton will 
travel through in a cylindrical system with a radius of 
a Pb nucleus. Such a scenario of weak interaction and 
small energy loss is barely consistent with the systemat- 
ics of the experimental data. It is clear that the large px 
spectra in Pb + Pb collisions at the CERN SPS energy 







can put very stringent limits on the interaction of ener- 
getic partons with dense medium and the induced energy 
loss. Within the parton model, one can exclude from the 
observed hadron spectra a parton energy loss larger than 
dE q /dx = 0.01 GeV/fm and a mean free path shorter 
than X q = 7 fm. This is much smaller than the most 
conservative estimate of parton energy loss in a dense 
medium [3-5]. According to the most conservative esti- 
mate in Ref. [5], a quark should still have an energy loss 
dE/dx 0.2 GeV/ fm in a cold nuclear matter of 10 fm in 
transverse size. Our constraint is 20 times smaller than 
this conservative limit, even though the dense matter in 
AA collisions should at least consists of hot hadronic mat- 
ter which is much denser than the normal cold nuclear 
matter. 

There are several implications one can draw from this 
analysis. Most of the recent theoretical estimates of par- 
ton energy loss are based on a scenario of a static and 
infinitely large dense parton gas. If the system produced 
in a central Pb + Pb collision only exists for a period 
of time shorter than the interaction mean free path of 
the propagating parton, one then should not expect to 
see any significant parton energy loss. Using the mea- 
sured transverse energy production dET/drj w 405 GeV 
[46] and a Bjorken scaling picture, one can indeed esti- 
mate [12] that the life time of the dense system in central 
Pb + Pb collisions is only about 2 - 3 fm/c before the 
density drops below a critical value of e c w 1 GeV/fm 3 . 
Even if we assume that a dense partonic system is formed 
in central Pb + Pb collisions, this optimistic estimate of 
the life time of the system could still be smaller than 
the mean free path of the propagating parton inside the 
medium. Thus, one does not have to expect a significant 
effect of parton energy loss on the final hadron spectra 
at large pr- Otherwise, it will be difficult to reconcile 
the absence of parton energy loss with the strong parton 
interaction which maintains a long-lived partonic system. 

Another conclusion one can also make is that the dense 
hadronic matter which has existed for a period of time 
in the final stage of heavy-ion collisions does not cause 
any apparent parton energy loss or jet quenching. A high 
Pt physical pion from jet fragmentation has a very long 
formation time. One does not have to worry about its 
scattering with other soft hadrons in the system which 
could cause suppression of high px pion spectra. We still 
do not understand the reason why a fragmenting parton 
does not loss much energy when it propagates through a 
dense hadronic matter. However, it might be related to 
the absence of energy loss to the quarks and anti-quarks 
prior to Drell-Yan hard processes in p + A and A + A col- 
lisions. This observation will make jet quenching a better 
probe of a long-lived partonic matter since one does not 
have to worry about the complications arising from the 
hadronic phase of the evolution. If one observes a dra- 
matic suppression of high pr hadron spectra at the BNL 
RHIC energy as predicted [2,7-9], then it will clearly in- 
dicate an initial condition very different from what has 
been reached at the CERN SPS energy. 



V. PREDICTIONS AT RHIC: IMPORTANCE OF 
P + A EXPERIMENTS 



As we have seen in the calculation of large pr particle 
spectra in p + p(p) collisions, the initial Ut smearing is 
most important at low energies where the original parti- 
cle spectra from jet production falls off very rapidly with 
Pt- The effect of fcy broadening due to multiple parton 
scattering in p + A collisions follows the same line. At 
higher energies {^/s > 50 GeV), when the original jet 
spectra become much flatter, the effects of both the ini- 
tial or intrinsic kx and nuclear kx broadening will become 
smaller, as we can clearly see in both our calculation and 
the experimental data in Figs. 8 and 5. At the BNL 
RHIC collider energy, the nuclear enhancement of large 
Pt spectra will then become smaller. Shown in Fig. 13 
are the calculated ratios of charged hadron spectra in 
p + Au over that in p + p normalized to the averaged 
number of binary nucleon collisions, 

p (n n_ d<j h AB /dyd 2 p T , . 

KabKPt) = 777 , , 1 „ — (Id) 

(N hinliI y)d<j£ p /dyd 2 pT 

where (iV binary ) = A for minimum-biased events of p + A 
collisions. We can see that the enhancement due to mul- 
tiple parton scattering will not disappear at the RHIC 
energy. There is still about 20-50% enhancement at mod- 
erately large pt around 4 GeV/c. It then disappears very 
quickly at larger pt . 

At the BNL RHIC energy, ^Jl = 200 GeV, nuclear 
modification of the gluon distribution will become im- 
portant for hadron spectra at moderately large pt- Such 
modifications for quark distributions have been measured 
in deeply inelastic lepton- nucleus collisions [47]. How- 
ever, the nuclear effects on the gluon distribution have 
not been directly measured. In the calculation shown as 
the dashed line in Fig. 13 we have used a recent parame- 
terization of the nuclear modification factors S u /a(x, Q 2 ) 
[Eq. (8)] by Eskola, Kolhinen and Salgado [47] which 
is based on global fits to the most recent collection of 
data available and some model on nuclear modification 
of gluon distribution. This result is quite different from 
the calculation using HIJING [38] parameterization. In 
EKS98 parameterization, QCD evolution equation has 
been used to take into account the Q 2 scale dependence 
of the nuclear modification which is absent in HIJING 
parameterization. Because of QCD evolution, the nu- 
clear shadowing effects at small x become smaller with 
increasing Q 2 . This is why the hadron spectrum in p + A 
collisions using EKS parameterization is larger than that 
using HIJING parameterization. The EKS98 parameter- 
ization also has a gluon anti-shadowing which is larger 
than any parameterizations before. Hadron spectra in 
a moderately large pt range at the BNL RHIC energy 
mainly come from fragmentation of gluon jets. This is 
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why the EKS98 result in Fig. 13 is larger than that with- 
out nuclear modification of parton distributions (solid 
line) in the px range where anti-shadowing becomes rel- 
evant. In addition, the impact parameter dependence of 
the nuclear modification of parton distributions is not 
implemented in the calculation using EKS98 parameter- 
ization. Though EKS98 parameterization has taken into 
account information about nuclear shadowing of gluon 
distribution from the measured scale evolution of the 
structure functions of different nuclei, such a procedure is 
also model dependent. The result will depend on whether 
one includes the higher gluon-density terms in the QCD 
evolution equation. Before a direct measurement of the 
nuclear effect on gluon distribution inside a nucleus, one 
should consider it one of the uncertainties in predicting 
the hadron spectra in p + A collisions at the BNL RHIC 
energy. 



1.6 



ravel any other effects in the hadron spectra caused by 
the formation of dense partonic matter in collisions. 
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FIG. 13. Predictions for the ratio of single charged hadron 
spectra in p + Au over p + p collisions normalized by the av- 
erage number of binary collisions (or A ) at i/s — 200 GeV. 
Different lines are for different parameterizations of shadow- 
ing or nuclear modification of parton distributions. 

As pointed out in Sec. II, the shape of the nuclear 
modification of hadron spectra at large pt is a result of 
multiple parton scattering inside a nucleus and their cor- 
responding absorptive corrections. The shape is remark- 
ably similar to the nuclear modification of parton distri- 
butions. In fact, nuclear modification of parton distribu- 
tions in small x region (shadowing and anti-shadowing) 
can be explained in a multiple scattering model [48]. 
Therefore, one could have double counted the same ef- 
fect if he includes both the Cronin effect (or kx broaden- 
ing) and nuclear modification of the parton distributions 
in the calculation of hadron spectra in p + A collisions. 
Given these uncertainties, it is therefore extremely im- 
portant to have a systematic study of p + A collisions at 
the BNL RHIC energy. Such an effort is pivotal to un- 
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FIG. 14. Predictions for the ratio of single charged hadron 
spectra in central Au + Au over p + p collisions normalized 
by the average number of binary collisions at ^/s = 200 GeV. 
Different lines are for different parameterizations of shadow- 
ing or nuclear modification of parton distributions. The upper 
set of lines are without parton energy loss and lower set are 
with parton energy loss dE q /dx = 0.25 GeV/fm and mean 
free path \ q — 2 fm. 



Shown in Fig. 14 are predictions for the ratio of hadron 
spectra in central Au + Au collisions over that in p + p 
[as defined in Eq. (16)], with and without parton energy 
loss. In both cases, uncertainties in the nuclear modifica- 
tion of parton distributions are still important. If there is 
significant parton energy loss in the dense medium at the 
BNL RHIC energy, leading hadrons from jets which are 
produced in the center of a overlapped region will be sup- 
pressed. As we can see this is where the Cronin effect and 
the effect of nuclear modification of parton distributions 
are the strongest. This is why the relative uncertainty 
caused by nuclear parton distributions is reduced com- 
pared to the case when there is no parton energy loss. In 
the calculations, we also included possible contribution 
to hadron production from soft processes at low pt as 
implemented in Ref. [9]. When hadrons from hard par- 
ton scattering are suppressed, these soft contributions 
might become dominant at intermediate px- This is why 
the spectra ratio increases again when px < 2 GeV/c in 
Fig. 14. The spectra at lower px could also be broadened 
due to the rescattering effect which also drives the sys- 
tem into local equilibrium. We should emphasize that the 
soft particle spectra we use here are extremely schematic 
and qualitative. This is another uncertainty one should 
keep in mind which could affect the shape of spectra at 
intermediate px- Nevertheless, the spectra at moderately 
large px is most sensitive to parton energy loss in dense 
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medium. Therefore, one should study hadron spectra in 
p + A very carefully and pin down the uncertainty due to 
multiple scattering and parton (anti)shadowing as much 
as possible. Only then one can draw more accurate con- 
clusions about parton energy loss from the single inclu- 
sive hadron spectra in heavy-ion collisions at the BNL 
RHIC energy. 



VI. CONCLUSIONS AND DISCUSSIONS 

In this paper, we have analyzed systematically large 
Pt hadron spectra in p + p, p + A and A + A collisions 
from CERN SPS to BNL RHIC energies within a QCD 
parton model. We found that both the initial kr in p + p 
collisions and the fcp broadening due to multiple parton 
scattering in p+A collisions are important to describe the 
experimental data within the parton model calculations. 
The value of initial kr in order to fit the data is found to 
be larger than the conventional value of 300 — 500 MeV 
for intrinsic kx according to the uncertainty principle. 
This finding is also consistent with analysis of Drell-Yan 
data [26] and recent study [20] of direct photon and pion 
production at around CERN SPS energy range. In both 
studies, one found that an intrinsic kr of the order of 1 
GeV/c is needed to describe the data within NLO par- 
ton model calculations. Since we only used LO pQCD 
calculation, we have to introduce some Q 2 dependence of 
the initial kr induced by initial-state radiation processes 
in order to fit the experimental data at low collider en- 
ergies. This parton model with phenomenologically in- 
cluded intrinsic kr describes very well the energy and 
isospin dependence of the hadron spectra. 

The parton model calculation can also describe the 
large pt pion spectra in heavy-ion collisions at the CERN 
SPS energies very well, both the A or centrality and en- 
ergy dependence. There is no evidence of proposed par- 
ton energy loss caused by dense partonic matter. Based 
on recent theoretical estimates [3-5] of parton energy loss 
in dense partonic matter, one should expect a parton en- 
ergy loss in the order of dE/dx ^2 — 4 GeV/fm. The 
absence of such energy loss in large pr hadron spectra 
implies that either there is no such dense partonic mat- 
ter formed or the life time of such medium is smaller 
then the mean free path of the parton interaction inside 
such a medium. It also tells us that the hadronic matter 
which must have existed for a period of time in heavy- 
ion collisions at the CERN SPS will not cause apparent 
energy loss or jet quenching effect. Therefore, if one ob- 
serves suppression of high pr hadrons at the BNL RHIC 
energy, it will unambiguously reflect an initial condition 
very different from what has been achieved at the CERN 
SPS. 

As pointed out in a recent paper by Gyulassy and Levai 
[49], even though HIJING [38] Monte Carlo model fails to 
reproduce the hadron spectra in p + p, p + A and S + S 
(and S + Au) data at the CERN SPS energy, it acci- 



dentally reproduces the central Pb + Pb data very well 
because some unique form of transverse momentum kick 
introduced in HIJING to the end-points of strings each 
time they suffer an interaction. A hydrodynamic model 
calculation [50] with significant transverse expansion can 
also be used to fit the high pt hadron spectra. So one 
should wonder how one can make sure that the high pt 
hadron production in A + B collisions is indeed domi- 
nated by hard parton scattering. As we have proposed 
earlier [12,51], measurement of two-particle correlation 
in azimuthal angle in the transverse plane should be able 
to distinguish these models from parton model. In the 
parton model, jets are always produced in pairs and back- 
to-back in the transverse plane. High pt particles from 
jet fragmentation should then have strong back-to-back 
correlation. Neither hydro-dynamical model nor multiple 
scattering of string end-points can give such correlation. 
If such correlation is seen, one can then study the pr de- 
pendence of the correlation to find out at what pr value 
the correlation disappears. One can then at least quan- 
tify above what pr value thermal-hydro model can be 
ruled out as the underlying particle production mecha- 
nism. Such a study is underway and will be reported 
elsewhere. 

At the BNL RHIC energy, things will become a lit- 
tle cleaner. As pointed out by Gyulassy and Levai [49], 
most of hadron production at moderate pt will be domi- 
nated by gluon production and are not influenced by the 
fate of strings or valence quarks (except net baryons). 
The spectra should be more sensitive to parton energy 
loss. Moreover, the hadron spectra in p + A and A + A 
collisions without parton energy loss will have a unique 
shape similar to those shown in Fig. 13 at high ener- 
gies which can not be fully revealed at the CERN SPS 
energies. This unique shape will be hard to explain by 
any thermal-hydrodynamic models. However, as we have 
demonstrated in the previous section, there are still some 
uncertainties related to nuclear modification of parton 
distributions. It is therefore very important to have a 
systematic study of p + p, p + A and A + A collisions 
at the BNL RHIC energy in order to make more quan- 
titative conclusions about parton energy loss in dense 
medium from the single hadron spectra at moderately 
high p T . 
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